Mon. Not. R. Astron. Soc. 000, [TH?] (0000) Printed 23 January 2013 (MN IATeK style file v2.2) 



Direct Gravitational Imaging of Intermediate Mass Black 
Holes in Extragalactic Halos 



o 

C 

(N 



o 

u 

6 

o3 



> 

O 
in 

o 



Kaiki Taro Inoue 1 * Valery Rashkov 2 f Joseph Silk 3 ' 4,5 J and Piero Madau 2 ; 

1 Department of Science and Engineering, Kinki University, Higashi-Osaka, 577-8502, Japan 

2 Department of Astronomy and Astrophysics, University of California, Santa Cruz, 1156 High Street, Santa Cruz, CA i 
3 Institut d'Astrophsique de Paris, UMR 7095, CNRS, UPMC Univ. Paris VI, 98 bis boulevard Arago, 75014 Paris, France 
4 Department of Physics and Astronomy, The Johns Hopkins University Homewood Campus, Baltimore, MD 21218, USA 

5 Beecroft Institute for Particle Astrophysics and Cosmology, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford, OX1 3RH, 



UK 



23 January 2013 



ABSTRACT 

A galaxy halo may contain a large number of intermediate mass black holes (IMBHs) 
with masses in the range of 10 2-6 M Q . We propose to directly detect these IMBHs by 
observing multiply imaged QSO-galaxy or galaxy-galaxy strong lens systems in the 
submillimctcr bands with high angular resolution. The silhouette of an IMBH in the 
lensing galaxy halo would appear as either a monopole-like or a dipole-like variation 
at the scale of the Einstein radius against the Einstein ring of the dust-emitting region 
surrounding the QSO. We use a particle tagging technique to dynamically populate a 
Milky Way-sized dark matter halo with black holes, and show that the surface mass 
density and number density of IMBHs have power-law dependences on the distance 
from the center of the host halo if smoothed on a scale of ~ 1 kpc. Most of the black 
holes orbiting close to the center are freely roaming as they have lost their dark matter 
hosts during infall due to tidal stripping. Next generation submillimeter telescopes 
with high angular resolution (< 0.3 mas) will be capable of directly mapping such 
off- nuclear freely roaming IMBHs with a mass of ~ 10 6 M Q in a lensing galaxy that 
harbours a O(10 9 ) M Q supermassive black hole in its nucleus. 
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1 INTRODUCTION 

Recent observations of off-nuclear ultraluminous X-ray 
sources (ULXs) suggest the presence of intermediate mass 
black holes (IMBHs) not only in the neighborhood of the 
galaxy nucleu s but also in star clu s ters far out in the 
galactic halo dMatsumoto et all 1200 ll : iRoberts et~afl 12004 

2010l ). A large population 



iFarrell et all 120091 : Ijonker et all 



of IMBHs might reside inside a galaxy halo, perhaps the 
leftover population of initial "seed" holes that never grew 
into the supermassive variety (SMBHfl) hosted today in 
the nuclei of massive galaxies. The mechanism of seed for- 
mation is unknown. Seed holes may be produced by the di- 
rect collapse of 10 4 — 10 6 Mq primordial gas clouds, by the 
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1 Here a "SMBH" means a black hole (BH) residing at the center 
of a main parent halo. If the host halo habouring the BH belongs 
to a more massive parent halo, we call it an "IMBH" . 



collapse of the first nuclear star clusters, o r be the rem- 
nants of 10 2 Mp) Population HI stars (e.g. | Loeb fc Rasiol 
1 1994 iMadau fc Reest l200ll : iDevecchi fc Volonteril 120091 ). If 
we could directly observe the abundance, spatial distribu- 
tion, and masses of IMBHs inside extragalactic halos, we 
would be able to constrain the process of SMBH seed for- 
mation that has hitherto been veiled in mystery. 

We propose here to directly detect IMBHs by observ- 
ing multiply imaged QSO-galaxy or galaxy-galaxy lens sys- 
tems in the submillimeter band with high angular resolu- 
tion (< 0.3 mas). Provided that the lensed source is suffi- 
ciently extended and the IMBHs inside the lensing galaxy 
are massive enough, then their gravita tional force would sig- 
nifica ntly distort the lensed images (|lnoue fc Chibal 1200a . 
l2006h . Moreover, any star cluster or dark matter subhalo 
that surrounds an IMBH would also distort the lensed im- 
ages at angular scales larger than the Einstein radius of 
the IMBHs. Using local distortions in the surface bright- 
ness of the lensed images, one can directly measure the po- 
sition, mass scale, and gravitational potential surrounding 
the IMBHs. 
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In this paper, we estimate the feasibility of mapping 
IMBHs in forthcoming observations. In Section 2, we de- 
scribe a particle tagging technique to dynamically populate 
the N-body Via Lactea II (VLII) extreme-resolution simu- 
lation with IMBHs, and derive semi-analytic formulae for 
describing the surface mass and number density of IMBHs 
in the host halo. In Section 3, we estimate the strong lens- 
ing probability due to IMBHs in a lensing galaxy. We then 
describe the extended source effects and observational fea- 
sibility of direct detection in the QSO-galaxy lensing sys- 
tem RXJ1131-1231. We summarise our results in Section 4. 
In what follows, we assume a concordant cosmology with a 
matter density Q. m = 0.272, a baryon density f2(, = 0.046, a 
cosmological constant S1a = 0.728, and a Hubble constant 
Ho — 70,km/s/Mpc, which are obtained from the observed 
CMB, the baryon acoustic oscillations, and measurement of 
H . 



2 A POPULATION OF IMBHS 

2.1 Simulation and tagging technique 

VLII, one of the highest-resolution N-body simulations of 
the assembly of a M ilky Way-sized galaxy halo to date 
IIDiemand et al ] |2008f ). provides an ideal set-up for the study 
of formation of IMBHs. VLII follows the hierarchical assem- 
bly of a dark halo of mass M200 = 2 x 10 12 Mq at redshift 
z = within r2oo = 402 kpc (the radius which encloses an 
average mass density 200 times the mean cosmological mat- 
ter density), with just over a billion particles and a force 
resolution of 40 pc. The simulat ion was perfo rmed with the 
parallel tree-code PKDGRAV2 (|Stadelll200lh . PKDGRAV2 
uses a fast multipole expansion technique in order to calcu- 
late the forces with hexadecapole precision, and an adaptive 
leapfrog integrator. Expected to harbor a thin disk galaxy, 
VLII was selected to have no major mergers after z = 1. 
At the present epoch, VLII has a cuspy density profile and 
exhibits rich galactic substructure - the main host's halo 
contains over 20,000 surviving subhalos with masses greater 
than 1O 6 M . 

Central black holes are added to subhalos following 
the pa rticle tagging technique detailed in lRashkov fc Madaul 
(|2013l ) and quickly summarized here. In each of 27 snapshots 
of the simulation, choosen to span the assembly history of 
the host between redshift z — 27.54 and the present, all 
subhalos are identified and linked from snapshot to snap- 
shot to their most massive progenitor: the subhalo tracks 
built in this way contain all the time-dependent structural 
information necessary for our study. We then identify the 
simulation snapshot in which each subhalo reaches its max- 
imum mass Mhaio before being accreted by the main host 
and tidally stripped. In each subhalo, the 1% most tightly 
bound dark matter particles are then "tagged" as stars at 
infall, following a set of simple prescriptions calibrated to 
reproduce the observed luminosity function of Milky Way 
satellites and the con centration of their stellar populations 
|Rashkov et a,l.|[2012l '). We then measure the stellar line-of- 
sight velocity dispersion, cr«, in each subhalo, and tag the 
most tightly bound central particle as a black hole of mass 
Mbh according to an extrapolation of the Mbh — cr* relation 



of lTremaine et alj (|2002r i. 



Mbh = 8.1 x 10 b M P 



(1) 

.100 km s" 1 / w 

For stellar systems with <r» ^ 6 km s _1 , a minimum seed 
hole mass of 100 M© is assumed. Any evolution of the tagged 
holes after infall is purely kinematical in character, as their 
satellite hosts are accreted and disrupted in an evolving 
Milky Way-sized halo. Black holes do not increase in mass 
after tagging, and are tracked down to the z — snapshot. 
We assume that IMBHs only form in subhalos with a mass 
at infall > M£ in = 10 7 M Q . 

The merging process that produces each subhalo prior 
to its infall into the main host halo is important, as black 
hole binaries may form in the process. The asymmetric emis- 
sion of gravitational waves produced during the coalescence 
of a binary black hole system is known to impart a velocity 
kick to the system that can displace the hole from the center 
of its host. The magnitude of the recoil will depend on the 
binary mass ratio and the direction and mag nitude of their 
spins, and follows the prescriptions detailed in lGuedes et alj 
(|201 ll ). When the kick velocity is larger than the escape 
speed of the host halo, a hole may be ejected into inter- 
galactic space before becoming a Galactic IMBH. IMBHs 
that would have formed in halos whose merger history points 
to a kick event are therefore e xcluded from the final catalog 
(see lRashkov fc Madaull2013l for details). 

Since some of the subhalos will fall on trajectories that 
bring them closer to the center of the main VLII host, they 
might be completely disrupted after infall, leaving freely- 
roaming ("naked") holes. Depending on the minimum halo 
mass that allows the formation of an IMBH, the numbers of 
formed, kicked, and stripped IMBHs will be different. Figure 
1 shows an image of the projected distribution of IMBHs in 
today's Galactic halo according to the above prescriptions. 
We find 1,070 naked IMBHs and 1,670 holes residing in dark 
matter satellites that survived tidal stripping. There are 50 
IMBHs more massive than 1O 3 M0, and about 2,500 holes 
at the minimum mass of 100 Mq (of which 950 are naked) . 
Naked holes are more concentrated towards the inner halo 
regions as a consequence of the tidal disruption of infalling 
satellites. Indeed, within 10 kpc, most MBHs are naked. 



2.2 Semi-analytic modeling 

To generalize the numerical results of the previous section, 
we shall assume in the following that the distribution of 
IMBHs inside a lensing galaxy halo is spherically symmet- 
ric, and that the surface mass density cr m and the surface 
number density a n are described by the following "universal 
profiles" , 



<7 m (r) 
(T n (r) 



O-m(O) 



((r/r c ) 2 + 1) exp[r/r,] ' 
M0) 



(2) 
(3) 



(r/r c + 1) exp[r/r»] ' 

where r c represents a core radius and r* denotes a cut-off ra- 
dius at which the density starts to decay exponentially with 
increasing r. The "universal profiles" are much steeper than 
those of a singular isothermal sphere (SIS). Subhalos hosting 
IMBHs can be massive enough to experience dynamical fric- 
tion, spiral in toward the center of the main host, be totally 
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Figure 1. Plot of the distribution of IMBHs in the simulated Milky Way-sized halo. Each dot represents the position of an IMBH (left). 
The threshold of the minimum subhalo mass required for the formation of a seed hole is assumed to be M*. = 10 7 Mq. The projected 
density of dark matter in the VLII simulation is plotted in a zoom-up image, where the size of each disk is proportional to the log 10 of 
the BH mass (right). 



stripped of their dark matter, and deposit a naked IMBH 
into the center of the main host. As shown in Figure 2, most 
of IMBHs at r < 10 kpc are naked. We find that the surface 
mass density at r < 10 kpc is not significantly affected by 
M s ■ 

JlJ mm- 

The constants o" m (0) and cr n (0) can be estimated as 
follows. First, each halo is approximated as a spherically 
symmetric object with virial radius r2oo- We assume that 
within r2oo the halo density profile is given by that of an 
SIS with one dimensional velocity dispersion a v . Then we 
have 



»"200(» 



(4) 



where Q m (z) and H(z) are the matter density parameter 
and the Hubble parameter at redshift z, respectively. The 
initial size of the halo is then 



200 



1/3 



T200- 



(5) 



The correlation between the stellar velocity dispersion a v 
and mass Msmbh of the supermassive black hole at the cen- 
ter is approximately given by 



/ Msmbh \ 



0x 



200 km/s 



(6) 



where /3 = Oil) and 4 < a < 5. The mass of the SMBH 
at the nucleus of the Milky Way is Msmbh = 4 x 10 6 M Q . 
Then, the velocity dispersion of the spheroidal component 
is a v = 88 km/s and r2oo = 3 .4 x 10 2 kpc at z = provided 
that a = 4.24 and /? = 1.32 (|Giiltekin et all 120091 ). We also 
assume that the total mass of IMBHs within 7-200 is given by 



Mimbh(< r 2 oo) = /Msmbh, where / = 0.1 - 1. Then from 
equations ([2]) and ((6J, we have 



.(0) 



Pf(a v /200 km/s) Q 
27rr 2 (ln(r 2 oo/r c ) + l/2] 



x 10 s M . 



(7) 



If we adopt / = 0.2, the analyticaly estimated surface mass 
density <r m fits the simulated values well (see Fig. 2). We 
find that the surface mass density in the neighborhood of 
the center does not change much even if one changes the 
threshold M^ in . 

Second, we assume that the seed of an IMBH is formed 
at a redshift of 2 = 20 inside a host halo with a mass of 
M > M,^ in , and each seed grows almost independently. Then 
the approximated virial radius of the halo is f 200(2 = 20) = 
3.3 x 10 2 pc. Assuming an SIS profile, the maximal circu- 
lar velocity of a host halo is estimated as Vmax = V2cr v = 
10//r2oo = 12 km/s. At z = 20, the number density of subha- 
los with maximal velocity larger than V ma x is approximately 
given by 



V 3 



A(z = 20) = 1.43 x 10 5 ( /i _1 Mpc/(km/s) 



(8) 



provided that Vmax < 1000 km/s (|Klypin et alj|201ll ). 

For the Milky Way-sized halo in VLII, the initial radius 
is estimated ELS Tint (2 = 20) = 1.65Mpc//i. Assuming that 
all the halos with mass > 10 7 Mq contain an IMBH seed, 
the total number of IMBH seeds inside a lensing galaxy halo 
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Figure 2. Surface mass density of IMBHs in the simulated Milky 
Way-sized halo. The dashed and dot-dashed curves correspond 
to the naked and the total (naked plus hosted-in-substructure) 
IMBHs in our numerical simulations, respectively. The solid 
curves represent fitted "universal profiles" in which the total 
mass of IMBHs coincides with that of our simulated results. This 
curve is also obtained if we assume / = 0.2. The distribution of 
simulated IMBHs is smoothed by a Gaussian window function 
W(r) = exp[— r 2 /(2rj?)], where r c = lkpc. The total masses of 
IMBHs (naked plus hostcd-in-substructure) are 9 X 1O 5 M0 (top) 
and 8 X 10 5 M o (bottom). 



is estimated as 

N(> V max = 12km/s) 



1.8 x 10 3 



d(£ = 2 °) x ^r 3 ■ 
12km/s) 3 3 



(9) 



Note that N(> Vmax) is proportional to rfoo- If the total 
number of IMBHs within r2oo does not change much, then 
the surface number density at the center is 



M0) 



N(> Vmax) 

27rr 2 oo?* c 



(10) 



Thus cr n (0) is proportional to r|oo or a% if r c does not de- 
pend on r2oo- As shown in Figure 3, the "universal profile" 
fits the simulation well if the total number of IMBHs coin- 
cides with that of the simulation. If a n (0) in equation (1101) is 
used, however, the surface number density is systematically 
reduced by a factor of 2 — 5 in comparison with the simula- 
tion. This difference may be due to the strong clustering of 
massive halos. 
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Figure 3. Surface number density of IMBHs in the simu- 
lated Milky Way-sized halo. The dashed and dot-dashed curves 
correspond to the naked and the total(nakcd plus hosted-in- 
substructure) IMBHs in our numerical simulations, respectively. 
The solid curves represent fitted "universal profiles" in which the 
total number of IMBHs coincides with that of our simulated re- 
sults (bold solid) and that of our semi-analytically estimated val- 
ues (thin solid). The distribution of simulated IMBHs is smoothed 
by a Gaussian window function W(r) = exp[— r 2 /{2rl)}, where 
r c = 1 kpc. 



3 LENSING EFFECT 
3.1 Lensing probability 

The Einstein angular radius 9e of a point mass with a mass 
M is written in terms of angular diameter distances: to the 
lens Dl, to the source Ds, and between the lens and the 
source Dls ss 



3 x 



M 



(l0 6 M o 



DlDs/Dl 
Gpc 



(11) 



Therefore, a radio interferometer with resolution of 3 mas 
can easily resolve the distortion of an image within the Ein- 
stein ring for a point mass M ~ 10 6 Mq. 

The strong lensing cross section due to an IMBH is pro- 
portional to M. Therefore, the lensing probability p is given 
by the ratio between the surface mass density of IMBHs 
and that of a lensing galactic halo at r — te — DlOe- From 
equations <|2j and (O, we have 



P : 



,(IMBH) 

m(SIS) 



r a — 4 

oc fa v , 



(12) 



where (T m (IMBH) and <r m (SIS) denote the surface mass den- 
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sity of IMBHs and that of an SIS, respectively. Similarly, the 
mean mass M of an IMBH at r = te satisfies 



M 



ff m (IMBH) 



ct„(IMBH) 



r a- 

oc fa v 



(13) 



Thus the lensing probability and the Einstein radius are 
larger for halos with larger velocity dispersion as long as 
a > 4 and equations © and ^ hold. As the Einstein radius 

— 1/2 

of an IMBH is proportional to D L , lens systems with 
small Dl are more suitable as targets. 



3.2 Extended source effect 

If a perturber is spatially extended, then the lensing effect 
is different from that of a point mass. The density profile 
of a perturber can be reconstructed from a local mapping 
between the observed image and the non-perturbed image 
obtained from the prediction of the macrolens. In fact, the 
power of the radial density profile of the perturbers can be 
reconstructed from the per turbed images within the Ein- 
stein ring of the perturber (jlnoue fc Chiball2005al ). In this 
way, one could make a distinction between an IMBH and a 
CDM subhalo. Furthermore, from distortion outside the Ein- 
stein ring of the perturber, the degeneracy between the per- 
turber mass and the distance can be broken provided that 
the Einstein radius of the perturbe r is sufficiently smalle r 
than that of the primary macrolens |lnoue fc Chiball2005bl ) . 
The precise measurement of spatial variation in the surface 
brightness of lensed images provides us with plenty of infor- 
mation about the mass, abundance, and spatial distribution 
of IMBHs. 



3.3 Simulation for RXJ1131-1231 

To estimate the observational feasibility of direct detection 
of IMBHs, we adopt a QSO-galaxy lensing system RXJ1131- 
1231 as a target since this system has a massive lensing halo 
at relatively small redshift. The redshifts of the source an d 
the lens are z s = 0.658 and z L = 0.295 (|Sluse et al.ll2003h . 
respectively. To model the macro-lens, we adopt a singular 
isothermal ellipsoid (SIE) in a constant external shear field 
in which the isopotential curves in the proje cted surface per- 
pend i cular to the line-of-sight are ellipses (|Kormann et al.l 
ll994l ; llnoue fc Takahashlll2012r ). The IMBHs inside the lens- 
ing halo are modeled as point masses. Using the observed 
mid-infrared fluxes, the position of the lensed images and 
the centroid of the lensing galaxy, the velocity dispersion 
for the best-fit model is estimated as a v — 3.5 x 10 2 km/s. 
From the M - a re lation with a = 4.24 and /3 = 1.32 
jGiiltekin et al.ll2009l '). the total mass of IMBHs within r 2 oo 
is Mimbh = 1-5/ x 10 9 M Q . Then the lensing probabil- 
ity is p = 5/ x 10~ 4 . The mean mass of the IMBH at 
r = r E = 5.6kpc/h is (Mimbh) = 2/ x 10 5 M Q and the 
corresponding Einstein radius is 8e = a/7 'mas. Therefore, if 
the angular resolution is < y^Jmas, wc would be able to de- 
tect an imprint of IMBHs in the lensed Einstein ring of dust 
emission if the lensed image has an area > 2000/~ 1 mas 2 . 

To estimate observational feasibility, we use the sim- 
ulated data of IMBHs for a Milky Way-sized halo with 
M^ in = 10 7 Mq and scale up the mass of each IMBH by 
a ratio between the total mass of the SMBH for RXJ1131- 
1231 (= 1.5 x 10 9 Mq) and that for the Milky Way-sized 



halo (= 4 x 1O 6 M0). In this model, we find that the ratio 
between the mass fraction of all IMBHs to that of a SMBH 
in the center is / — 0.2. The most massive IMBH has a mass 
of Mimbh = 7 x 10 r Mq . Taking into account the logarith- 
mic correction to <j v in equation (0, the distance of each 
particle from the center is scaled up by a factor 



ln(a„(RXJ1131)/(V50ff(,z = 0.295)r c )) 



ln(<7„(MW)/(V50//(z = 0)r c )) 



1/2 



= 1.05, 



(14) 



where we assume that r c = 1 kpc. As shown in Figure 4, 
we find that massive "naked" IMBHs with masses in the 
range of O(1O 5 " 6 )M are observable if the radius of the 
dust emitting region around the QSO is ~ 500 Mpc and the 
angular resolution is < 1 mas. The area of the lensed arc is 
found to be ~ 6x 10 5 mas 2 . Thus we expect > O(10) "naked" 
IMBHs inside the lensed arc if / > 0.2. 

Even if IMBHs reside in places far from the brightest 
part of the Einstein ring of dust emission, the surface bright- 
ness of the lensed image that is sufficiently close to the posi- 
tion of an IMBH is the same as that of the brightest region 
in the Einstein ring. It is easily distinguished from other 
sources as the spatial gradient of the surface brightness is 
extremely large near the position of an IMBH. 



4 SUMMARY 

Based on numerical simulations of IMBH formation in a 
Milky-way sized halo, we have found that the surface mass 
density and number density of IMBHs have power-law de- 
pendences on the distance from the center of the host halo 
if smoothed on scales of ~ 1 kpc. Most of the black holes 
orbiting close to the center are freely roaming as they have 
lost their surrounding dark mass during infall due to tidal 
stripping. Assuming the surface mass density and number 
density have such power-law dependences, the strong lens- 
ing probability due to free roaming IMBHs in a massive 
lensing galaxy such as RXJ1131-1231 is ~ O(10~ 4 ). The 
next generation submillimeter telescopes with high angular 
resolution (< 0.3 mas) will be capable of directly mapping 
these freely roaming IMBHs with a mass of ~ 10 6 M . 

In addition to IMBHs in the lensing galaxy, SMBHs 
with a mass of > 10 6 Mq inside halos in the line-of-sight may 
also be observable, especially in systems that show anomalies 
in the flux ratios (|lnoue fc Takahashll2012l ). In this case, 
some distortion in the lensed image due to the host halo is 
expected. 

By measuring the local distortion of lensed extended 
images with high angular resolution, we will be able to deter- 
mine the mass, abundance, and density profile of the IMBHs 
present in the lensing galaxy. Direct detection of IMBHs 
will shed new light on the formation process of SMBH seeds 
which has hitherto been shrouded from view. 
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Figure 4. Simulated lensed images of RXJ1131-1231. The black dots in the top left panel show the positions of IMBHs in the lensing 
galaxy and the dot size is proportional to log 1 g(Afi]viBH)- The other panels show the contour and 3D maps of the surface brightness 
centered at the angular position of an IMBH with a mass of 1.4 X 10 6 Mq that produces a "black hole". We assume that the dust 
emitting region has a circular gaussian luminosity profile with a 1 a radius r = 500 pc. The numbers of contours representing the surface 
brightness are 10 (top left), 25 (bottom left), 25 (top right). 



(DARK) hosted by Universite Pierre et Marie Curie - Paris 
6. 



REFERENCES 

Devecchi B., Volonteri M., 2009, The Astrophysical Jour- 
nal, 694, 302 

Diemand J., Kuhlen M., Madau P., Zemp M., Moore B., 

Potter D., Stadel J., 2008, Nature, 454, 735 
Farrell S. A., Webb N. A., Barret D., Godet O., Rodrigues 

J. M., 2009, Nature, 460, 73 
Guedes J., Madau P., Mayer L., Callegari S., 2011, The 

Astrophysical Journal, 729, 125 
Giiltekin K., Richstone D. O., Gebhardt K., Lauer T. R., 

Tremaine S., Aller M. C, Bender R., Dressier A., Faber 



S. M., Filippenko A. V., Green R., Ho L. C, Kormendy 
J., Magorrian J., Pinkney J., Siopis C, 2009, The Astro- 
physical Journal, 698, 198 

Inoue K. T., Chiba M., 2003, The Astrophysical Journal 
Letter, 591, L83 

Inoue K. T., Chiba M., 2005a, The Astrophysical Journal, 
634, 77 

Inoue K. T., Chiba M., 2005b, The Astrophysical Journal, 
633, 23 

Inoue K. T., Chiba M., 2006, Annual reports by Research 
Institute for Science and Technology, 18, 11 

Inoue K. T., Takahashi R., 2012, Monthly Notices of Royal 
Astronomical Society, 426, 2978 

Jonker P. C, Torres M. A. P., Fabian A. C, Heida M., 
Miniutti C, Pooley D., 2010, Monthly Notices of Royal 



Direct Gravitational Imaging of Intermediate Mass Black Holes in Extragalactic Halos 7 



Astronomical Society, 407, 645 

Klypin A. A., Trujillo-Gomez S., Primack J., 2011, The 
Astrophysical Journal, 740, 102 

Kormann R., Schneider P., Bartelmann M., 1994, Astron- 
omy and Astrophysics, 284, 285 

Loeb A., Rasio F. A., 1994, The Astrophysical Journal, 
432, 52 

Madau P., Rees M. J., 2001, The Astrophysical Journal 

Letter, 551, L27 
Matsumoto H., Tsuru T. G., Koyama K., Awaki H., 

Canizares C. R., Kawai N., Matsushita S., Kawabe R., 

2001, The Astrophysical Journal Letter, 547, L25 
Rashkov V., Madau P., 2013, in preparation 
Rashkov V., Madau P., Kuhlen M., Diemand J., 2012, The 

Astrophysical Journal, 745, 142 
Roberts T. P., Warwick R. S., Ward M. J., Goad M. R., 

2004, Monthly Notices of Royal Astronomical Society, 349, 

1193 

Sluse D., Surdej J., Claeskens J. F., Hutsemekers D., Jean 
C, Courbin F., Nakos T., Billeres M., Khmil S. V., 2003, 
Astronomy and Astrophysics, 406, L43 
Stadel J. G., 2001, Ph.D. Thesis, University of Washington 
Tremaine S., Gebhardt K., Bender R., Bower G., Dressier 
A., Faber S. M., Filippenko A. V., Green R., Grillmair 
C, Ho L. C, Kormendy J., Lauer T. R., Magorrian J., 
Pinkney J., Richstone D., 2002, The Astrophysical Jour- 
nal, 574, 740 



